نشریه مهندسی عمران فردوسی سال بیست و پنجم شماره دی ۱۳۹۳ 


ارائه‌ی الگور یتم‌های بهینه‌سازی کارا با استفاده از روش رشد سا ز کار پیولوژیکی * 
هی وی ی تس ری 


چکیده تاکنون روش‌های محتلفی برای بهینه‌سازی شکل و ساتار سازه‌های مکانیکی با الهام از رفتار طبیعت تدوی نگردیده است. در این 
مقاله با استفاده از اصل رشد سازگار بیولوژیکی, دو الکوریتم برای بهینه‌سازی شکل و ساختار ارائه شده است. الگوریتم اول به یک صفحه‌ی 
سوراخ دار و نیز یک فیلت و الکوریتم دوم به یک تیر مخروطی و نیز یک تیر با سطح فطع متفاوت اعمال شده است. نتایج نشان می‌دهد که 
عمر خستگی صفحه‌ی سوراخ دار تا ۶۰ برابر افزایش و تنش‌های ماکزیمم فیلت نسبت به تتایج دیگر پژوهشگران "درصد کاهش يافته است . 


هم‌چنین این الکوریتم می‌تواند حجم سازه را به طور قابل توجهی کاهش دهد. 


واژه‌های کلیدی بهینه‌سازی شکل. بهینه‌سازی ساختان رشد سازگار بیولوژیکی. عمر خستگی. 
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تاریخ دریافت مقاله ٩۰/۹/۹‏ و تاریخ پذیرش آن ۹۲/۷/۱۳ می‌باشد. 
(۱) کارشناس ارشد مکانیک دانشگاه فردوسی مشهد. 
(۲) نویسنده‌ی مسول استاد گروه مکانیک. مرکز پژوهشی مهندسی تولید ناب دانشگاه فردوسی مشهد. 


‌ 


معد مه 


بهینه‌سازی شکل و ساختار سازه‌های مکانیکی یکی از 
موضوعاتی است که در سال‌های اخیر مورد توجه 
بسیاری از پژوهشگران بوده و روش‌های مختلفی برای 
آن ازاته شده است: تروش‌هاق بهینه‌سازی شکل سبولا 
یک ناحیه‌ی کوچک از جسم را مورد مطالعه قرار 
می‌دهد و عملیات بهینه‌سازی تنهاروی این قسمت 
انجام می‌گیرد تا برای مثال ماکزیمم تنش موضعی در 
آن قسمت مینیمم شود.در اين نوع از بهینه‌سازی هنگام 
شروع عملیات ساختار جسم به طور کامل مشخص 
است و قسمت‌هایی از سازه مثلاً سوراخ‌ها و یا 
شکاف‌ها دارای تنش‌های موضعی غیر یکنواخت 
می‌باشند. هدف این روش‌ها پیدا کردن بهترین شکل و 
پروفیل برای قسمت‌های ذکر شده است تا ماکزیمم 
تنش به اندازه‌ی کافی کوچک شود و مرز مورد مطالعه 
دارای تتش یکنواخت‌تری گردد. 

سه روش اصلی برای بهینه‌سازی شکل وجود 
دارد. این روش‌ها عبارتند از: برنامه‌نویسی ریاضی [1] 
روش رشد سازگار بیولوژیکی[2]و بهینه‌سازی تکاملی 
سازه‌ها [3]. روش برنامه‌نویسی ریاضی بر پایه‌ی 
ریاضیات استوار است و برای محاسبه‌ی شکل نهایی 
بهینه به محاسبه‌ی گرادیان تنش نیاز است. برای 
محاسبه‌ی گرادیان تنش ابتدا گرادیان تغییر مکان نسبت 
به متغیرهای طراحی محاسبه می‌شود و سپس با استفاده 
زو اک زب رتخا کزامتان کش فتاسته 
می‌گردد. پس از محاسبه‌ی گرادیان‌ها. برای محاسبه‌ی 
متغیر طراحی می‌توان از هر روش ریاضی غیر خطی 
برای حل معادلات به دست آمده استفاده کرد. 
بهینه‌سازی شکل با استفاده از برنامه‌نویسی ریاضی 
روش بسیار قدرتمندی است ولی دارای اشکالات 
میاستن لت اتتر کترین افسکال اه روش تم 
محاسبات بالا و در نتیجه هزینه‌ی بالای این روش 
ست. از طرفی به علت انجام عملیات ریاضی بی‌شمار 
برای محاسبه‌ی گرادیان‌ها و متغیرهای طراحی: دقت 
ین روش در قیاس با روش‌های دیگر مانند روش رشد 


نشریه مهندسی عمران فردوسی 


ارائه الگوریتم‌های بهینه‌سازی کارا با استفاده از روش رشد ... 


سازگار پیولوژیکی کم‌تر است [4], 

از مسا رم هار زای استفاهه هر کته 
هیچ محدودیتی در شکل نهایی سازه یا قطعه‌ی 
مکانیکی وجود ندارد و شکل قطعه تا پایان فرایند 
بهینه‌سازی نامشخص است. تنها اطلاعات موجود. 
معیار بهینه‌سازی و محدودیت‌های طراحی از قبیل 
نیروها و تکیه‌گاه‌ها می‌باشد.هدف بهینه‌سازی ساختار 
حذف يا باز توزیع ماده در دامنه‌ی اولیه‌ی سازه و 
رسیدن به یک ساختار بهینه با توزیع تنش نسبتاً 
یکنواخت تری می‌باشد. این روش یک روش مبتنی بر 
تکرار می‌باشد و معمولاً همراه با تاریخچه‌ای از 
ساختارهای میانی است که آن‌ها نیز در فرآیند طراحی 
قابل استفاده‌اند. 

تاکنون چهار الگوریتم اصلی برای بهینه‌سازی 
ساختار سازه‌ها ارائه شده است [5]. ایین الگوریتم‌ها 
عبارتند از: روش حسذف سخت الرم. روش 
همگن‌سازی. دیدگاه سازه مبنا و روش حباب. در سه 
روش آخر باید تابع هدف. متغیرهای طراحی و قیود را 
تعریف نمود مسئله‌ی بهینه‌سازی را با استفاده از 
روش‌هایی مانند برنامه‌ریزی ترتیبی درجه‌ی دو حل 
کرد. این سه روش نسبتاً پیچیده‌اند.ولی بر خلاف آن‌ها 
روش حذف سخت انرم بر اصول بسیار ساده‌ای استوار 
است. در این پژوهش نیز از الگوریتم حذف نرم ( 5011 
(8160) عمنام0 النکاستفاده شده است. 


روش رشد سازگار بیولوژیکی در 
بهینه‌سازی شکل 
نظریه‌ی این روش بر اساس تغییرات سازه‌های 
بیولوژیکی در اثر نیروهای مختلف در طبیعت می‌باشد. 
در این روش نبازی به استفاده از گرادیان‌های تنش و 
کت عیرست تن اسان سر اه 
محاسبه‌ی گرادیان باعث می‌شود که سرعت عملیات 
بهینه‌سازی بسیار زیاد و هزینه‌ی بهینه‌سازی به طور 
قابل توجهی کاهش یابد. دقت این روش بسیار بیش‌تر 
از سایر روش‌های بهینه‌سازی شکل است. به گونه‌ای 
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که بعضی از پژوهشگران از نتایج حاصل از روش رشد 
سازگار بیولوژیکی به عنوان معیاری برای ارزیابی 
روش‌های دیگر در تحقیقات خود استفاده کرده‌اند 
[6-8]. نظریهی بنیادی در روش رد سازگار 
بیولوژیکی این است که برای مینیمم سازی تنش‌های 
ماکزیمم در محل‌هایی که تنش موضعی زیادی وجود 
دارد ماده اضافه گردد و در محل‌هایی که تنش موضعی 
کم‌تر است از مقدار ماده کاسته شود. این عملیات با 
تغییر موقعیت گره‌ها در مدل المان محدود انجام 
می‌شود و تا جایی ادامه می‌یابد که مقدار تتش در مرز 
مورد مطالعه یکنوانحت شود و شرط همگرایی برآورده 
گردد. 

برای بهینه‌سازی به روش فوق پژوهشگران 
الگوریتم‌های مختلفی ارائه کرده‌اند[9-11]. الگوریتمی 
که در این پژوهش مورد بررسی قرار گرفته الگوریتم 
پیشنهادی [۷۷۷]2 است. عملکرد این الگوریتم که بر 
پایه‌ی اختلاف تنش‌های ون مایزز بنا گردیده. بدین 
صورت است که برای به دست آوردن شکل بهینه‌ی 
یک سوراخ باید اختلاف تنش‌های ماکزیمم و مینیمم 
وزساپوز از یک تیار اس کر توافت ان 
لکرتخ در تن ی ریم بر شود کل زز بر 


مراحل ف رآیند بهینه‌سازی.در این پزوهش الگوریتم 

پیشنهادی [۷]2به صورت اصلاح شده‌ی زیر به‌کار 

گرفته می‌شود: 

۱- تعریف پارامترهای طراحی و بهینه‌سازی شامل: نوع 
المان به‌کار گرفته شده برای تحلیل اجزای محدود. 
خواص مواد. ابعاد. بارگذاری‌ها؛ شرایط مرزی؛ 
فاکتور مقیاس, تعداد نقاط کنترلی؛ تعداد تکرار و 
دفت بهینه‌سازی. 

۲- مدل‌سازی اولیه و تحلیل اجزای محدود مدل اولیه: 
این مرحلسه با استتفاده از نسرم‌افزار 
انسیس(۸(۹۷۹) [12]انجام می‌گیرد. 


سال بیست و پنجم شماره دو. ۱۳۹۳ 


۳- هم گرایی بهینه‌سازی: در این مرحله پارامتر تلبه 
صورت معادله‌ی (۱) تعریف می‌شود: 


0 


که سس و ,م6 به ترتیب تنش‌های ون مایزز ماکزیمم 

و مینیمم هستند. حال اگر مقدار از دقت تعیین شده 

برای فرایند بهینه‌سازی (7) کم‌تر باشد. عملیات متوقف 

می‌گردد. در غیر این صورت وارد چرخه می‌شود. 

۶- چرخه : در این مرحله باید مختصات جدید نقاط 
کنترلی را از معادله‌ی (۲) به دست آورد. 


۳( (ز60(<1 +( << (260 


که در آن (26)1 و ([)به ترتیب مختصات جدید و قبلی 
نقطه‌ی ام است. () از معادله‌ی (۳) برای نقطه 
کنترلی 1 به دست می‌آید و (ز)16فاکتور مقیاس در 
چرخه‌ی ژام است و مقدار اولیه‌ی ۰/۰۰۰۱ برای آن در 


نظر گرفته می‌شود. 
۳ (رر6/((ف) 6 - ()بری6) ع (فاه 


که () یم نش ون مایزز در نقطه‌ی کترلی ذام است, 

۵- مدل‌سازی مجدد و تحلیل اجزای محدود مدل 
جدید: در این مرحله سازه بار دیگر با استفاده از 
نقاط کنترلی جدید مدل‌سازی می‌شود. سپس 
تحلیل المان محدود انجام می‌گردد. 

1- هم گرایی بهینه‌سازی: در این مرحله اگر ۲۳ در 
چرخه‌ی رام کم‌تر از چرخه‌ی 1-زام باشد» 
فاکتور مقیاس در عدد ۰/۵ و در غیر این صورت 
در عدد ۱/۱ ضرب می‌شود. 

۷- مراحل ۳ تا 1 تا آن‌جا تکرار می‌شود که مقدار :8 


از دقت بهینه‌سازی (7) کم‌تر شود. [2] 


نشریه مهندسی عمرآن فردوسی 


۲۰ ارائه لگوریتم‌های بهینه‌سازی کارا با استفاده از روش رشد ... 


ساخت مدل هندسی 
تعیین مسیر حرکت نقاط انتهایی 


شکل ۱ الگوریتم پیشنهادی [۷۷1]2برای فرایند بهینه‌سازی شکل 


100 ۵ 100 ۵ 


1 ۶ 


2 


۱0۱۱۱۹۱۱ 


شکل ۲ صفحه دارای سوراخ مرکزی الف) با بارگذاری کششی. ب) با بارگذاری فشاری 
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کاربرد روش رشد سازگار بیولوژیکی در 
بهینه‌سازی شکل 

مثال ۱-یک صفحه‌ی مربعی دوبعدی دارای سوراخ 
مرکزی که از دو ظرف تحت بارگذاری نامتقارن قرار 
گرفته در نظر بگیرید. هدف بهینه‌سازی شکل صفحه به 
منظور افرايش عمر خستگی آن است (شکل ۲- الف یا 
ب). طول ضلع مربع. قطر سوراخ و ضخامت صفحه به 
ترتیپ 2۰۰ 1*۰ و *عمیلی‌متر در نظر گرفته شده‌اند. 
همم‌چنین در تسام تحلیل‌هاماده با مدول 
الاستیسیته‌ی ۲۱۰۳۵ و ضریب یواسون ۰/۳ فرض شده 
ای گرا اال بان ری ات متس و سار 
کششی در جهت / به میزان ۱۰۰۸۲۳۵ و در جهت ۶ به 
۵۱/۳۵ سک (2۷ الف او سس کته بیان 
فشاری به همین میزان شکل (۲- ب) قرار می‌گیرد. 
فرک‌انس بار نوسانی برابر 107016/5در نظر گرفته 
ی و3 

به دلیل تقارن. یک چهارم صفحه در نرم‌افزار 
5 تحلیل شده است. نتایج شکل (۳) نشان 
می‌دهد که در اطراف سوراخ به دلیل تمرکز تنش و 
توزیم غیریکنواخت آن ترک ناشی از خستگی به 
سرعت رشد می‌کند و عمر سازه در اطراف سوراخ به 
شدت کاهش می‌یابد. 


شکل ۳ توزیع غیریکنواخت تنش در اطراف سوراخ بهینه نشده 
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۳۱ 


ق رآیند بهینه‌سازی. بر مبنای الگوریتم بیان شده با 
استفاده از زبان طراحی پارامتری انسیس (۵- ۸9۷5 
(اناطه) معمتاعصها معلوهنا متتامصصهت۳2) کدی 
نوشته شده که شکل بهینه‌ی سوراخ را مشخص می‌کند. 
همان‌گونه که در شکل (۶) مشاهده می‌شو شکل 
بهینه شبیه یک بیضی است. به گونه‌ای که قطر بزرگ 
بیضی در امتداد نیروی بزرگ‌تر و قطر کوچک در 
امتداد تیروی کرهکقر استهو ایس با ال رش 
داز کان زوین که زرف ایا کیش وان 
ماده افزوده شود و در مقاطع با تنش کم‌تر از مقدار ماده 
کاسته شود. کاملا مطابقت دارد. 


شکل؟ توزیع یکنواخت تنش در اطراف سوراخ بهینه 


در شکل (۵) برای نمایش واضح‌تر, اطراف سوراخ 
با بزرگ‌نمایی نشان داده شده است. همان‌گونه که در 
شکل‌های (۵و 4) دیده می‌شود. توزیع تنش در اطراف 
سوراخ یکنواخت است و تش‌ها در محدوده‌ی 
کوک فری گرا دز مان ای سا کزیی اما 
یافته و تمرکز تنش بسیار کم‌تر است.پارامتر 158 به 
عنوان درصد کاهش تنش در حالت بهینه نسبت به 
حالت اولیه به صورت معادله‌ی(4) تعریف می‌شود. 


135 < ])6, - (0 (2 


که :6و ,6 به ترتیب تنش‌های اولیه (غیر بهینه) و بهینه 


نشریه مهندسی عمرآن فردوسی 


۳۲ 


می‌باشند. اگر در معادله‌ی (4) به جای :0و ,0 به 
ترتیب تلاو ,1 قرار داده شود درصد افزایش 
یک‌نواختی تنش به دست می‌آید که در آن :18۳و ,13۳ 
از معادله‌ی (۱) به ترتیب برای حالات بهینه نشده و 
بهینه به دست می‌آیند. در جدول (۱)تنش‌های ون 
مایزز در دو حالت بهینه نشده و بهینه مقایسه شده 
است. نتایج نشان می‌دهد که تنش ون مایزز متوسط 
بترا ات تشه ۲۲ فرش کاهق یافته لس تا: 
هم‌چنین مقدار 18 برای 47,187 درصد به دست می‌آید 
که نشان می‌دهد تنش در حالت بهینه ۶۷ درصد 


یک‌نواخت‌تر شده است. 


جدول اتنش‌های ون مایزز در دو حالت بهینه نشده و بهینه برای 


صفحه‌ی دارای سوراخ مرکزی 


۳9 (۸۲۰) ۹5 ۷۵۳۰965 8 
صیاصمته]۱۷۱ ز ‏ مصیاوصت از صهع]۷ 

0.607 | 1.50422 2 2.262 ۳ 

نشده 

بهینه 133912 2 ۳ 02 0.36 

13 

47 22 37 0۳0 


در نمودار شکل (1) نتایج این پژوهش با نتایج 
|۱2 یه قنله انیت پراش ماس قق ها رون 
مایزز ماکزیمم نرمال شده مربوط به هر چرخه از 
معادله‌ی (۵) استفاده می‌شود. 
۹ 0۵(1 - مس - مسد سه 
که سرد سس مه به ترتیب تنش ون مایزز ماکزیمم 
و نرمال شده در چرخه‌ی [ ام و ,,6 تنش مرجع است. 
م0 را می‌توان تنش ون مایزز صفحه بدون سوراخ در 
نظر گرفت که برایر ۸۷/۱۸۳۵ است. در نمودار شکل 
(1) سرعت همگرا شدن الگوریتم به‌کار برده شده 
مشاهده می‌شود. با مشاهده‌ی نمودار می‌توان نتیجه 
گرفت که بعد از ۲۰ چرخه جواب همگرا شده است و 
این با نتایج [2]/کاملاً سازگار است. همان‌گونه که از 


نشریه مهندسی عمران فردوسی 


ارائه الگوریتم‌های بهینه‌سازی کارا با استفاده از روش رشد ... 


نمودارها مشخص است. تنش‌های ماکزیمم نرمال شده 


تقریباً با هم برابرند. 


شکل ۵ توزیع یکنواخت تنش در اطراف سوراخ بهینه با 


بزرگ‌نمایی اطراف سوراخ 


ِ‌ 

1 
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شکل 1 سرعت همگرا شدن الگوریتم به کار پرده شده برای 
صفحه‌ی سوراخ دار الف) در این پژوهش. 
ب) الگوریتم [۲۷۷۷]2 


محاسیه‌ی عمر خستگی.در اثر اعمال بار تناوبی بر 
روی سازه. با این که تنش ماکزیمم ایجاد شده کم‌تر از 
تنش نهایی آن است. اما پس از تعدادی سیکل, بر روی 
سازه ترک‌هایی ایجاد شده که در نهایت منجر به 
شکست می‌شود. به این پدیده خستگی در اثر اعمال 
بار تناوبی. و به تعداد سیکل‌هایی که سازه قبل از 
شکست تحمل می‌کند. عمر خستگی می‌گویند. روش 
محاسبه‌ی عمر خستگی در نرم‌افزار انسیس بر اساس 
استاندارد مخازن تحت فشار (۷۵55615 ۳۳۵55۱۲6) 


است. در این پژوهش هدف بهینه‌سازی شکل صفحه 


سال بیست و پیجم» شماره دی ۱۳۹۳ 


ترا اقزایتن عم عسعگی ان استه بسن از بسهخست 
آمدن شکل بهینه‌ی صفحه‌ی دارای سوراخ مرکزی باید 
عمر خستگی در اثر بارهای تناوبی کششی و فشاری 
برای تمامی گره‌های مدل محاسبه گردد. برای 
محاسبه‌ی عمر خستگی توسط نرم‌افزار انسیس باید 
گام‌های زیر برداشته شود: 
۱- مدل‌سازی و حل مسئله با توجه به بارگذاری 
موجود در حل خحستگی: قبل از انجام هر تحلیل 
شستگی باید یک تطیل استاتیکی که.شامل حداقل دو 
بارگذاری است انجام شود؛ که در اینن پزوهش 
صفحه‌ی با دو بارگذاری کششی و فشاری تحلیل 
استاتیکی می‌شود. 
۲- تعیین خواص خستگی ماده‌ی مورد استفاده شامل 
نمودار تتش سیکل (0۷70169 01 06۲ _ووع0٩‏ ). 
۳- تعیین خواص تست کشش- فشار ماده مورد 
استفاده شامل تتش نهایی (5۷» کاهش سطح 3۸ 
مدول الاستیسیته (۲). 

در اینجا یز با استفاده از مآللهکدی نوشته شسده 
تا بتواند عمر خستگی صفحه را به صورت توزیع عمر 
برای تمامی گره‌های مدل بر روی مدل نمایش دهد. در 
ابتدا این کد به صفحه‌ی بهینه نشده اعمال می‌شود. 
نتایج عمر خستگی به صورت شکل (۷) بر روی 
صفحه نشان داده شده است. همان گونه که مشاهده 
می‌گردد در گوشه‌ی پایین سوراخ عمر صفحه برابر 
۸ سیکل است. بنابراین صفحه پس از ۵۰۸ سیکل 
کشش و فشار از گوشه‌ی پایین سوراخ ترک برمی‌دارد 
و شروع به شکست می‌نماید. 

حال اگر کد نوشته شده به صفحه با سوراخ بهینه 
شده اعمال گردقه مشاهده می‌شود که عمر خستگی 
بسیار افزايش می‌یابد. نتایج عمر خستگی به صورت 
شکل (۸) بر روی صفحه نشان داده شده است. 
همان‌گونه که در شکل دیده می‌شود. عمر صفحه برابر 
۱ سیکل است. یعنی در حالت بهینه. صفحه‌ی 
دیگر از گوشه پایین سوراخ ترک برنمی‌دارد و پس از 
۷ سیکل کشش و فشار به تدریج از گوشه‌ی 


سال بیست و پنجم. شماره دو ۱۳۹۳ 


۳۳ 


بالای سوراخ بیضوی ترک بر می‌دارد. 

مقایسه‌ی بین دو حالت نشان می‌دهد که در حالت 
پهینه‌ی عمر صفحه حداقل ۶۰ برابر افزايش می‌یابد و 
این در حالی است که هنوز آثار ترک چندان به چشم 
نمی‌آید. مسئله‌ی مورد بررسی این مقاله (محاسبه‌ی 
عمر خستگی) و روش بهینه‌سازی استفاده شده در 
مراجع دیگر مورد بحث قرار نگرفته است. با اين همه 
در پژوهش متک و دیگران [13] نیز عمر خستگی 
سازه‌ه ا با روش بهینه‌سازی به وسیله‌ی کامپیوتر 
((0۸0) عماح‌عتصصتامون 060همانام‌ططهن)) تا *ع 


شکل ۷ نتایج عمر خستگی برای صفحه‌ی بهینه نشده و کاهش 


مقدار آن در گوشه‌ی پایین سوراخ 


شکل ۸ نتایج عمر خستگی برای صفحه‌ی بهینه و افزایش آن 
نسبت به حالت بهینه نشده 


مثال ۲- در این مثال یک فیلت دوبعدی که از یک 


طرف تحت تنش یک‌نواخت و از طرف دیگر روی 


نشریه مهندسی عمران فردوسی 


۳ 


تکیه گاه ساده قرار گرفته است؛ برای یک نواخت سازی 
توزیع تنش در محل ناچ بررسی می‌شود(شکل4). این 
مسئله نیز توسط [۷۷]2آبررسی شده است. مشخصات 
سازه عبارت اسشحرن از: 0۳۵ 210 < ۳ 0.3 < ۰10 
0.5 < 5 < 51 < با بل 03 < ۵ 5 02 < 9 
کته عر ان ۲ مدول الاستیسیته و 9 ضریب پواسون 
است. به علت تقارن تنها نیمی از فیلت با استفاده از 
نرم‌افزار 5 تحلیل شده که در شکل (۱۰) نشان 
داده شده است. 

برای اعمال الگوریتم رشد بیولوژیکی با استفاده از 
, کدی نوشته شده که شکل بهینه‌ی فیلت را بعد 
از ۲ مرتبه تکرار به دست شا روگ توزیع تنش ود 
مایزز برای حالت‌های بهینه و غیر بهینه در شکل‌های 
(۱۲و۱۱) ارائه شده است. همان‌گونه که مشاهده 
می‌شود در حالت بهینه شعاع فیلت برای از بین بردن 
تنش‌ها در بازه‌ی کوجک‌تری توزیع شده‌اند. نتایج 
نشان می‌دهد که مقدار 16 برای تنش ون مایزز 
متوسط ۰ درصد. اسینتا: در حدول (۲) تنش‌های ون 


مایزز در حالت‌های بهینه و غیر بهینه مقایسه شده 


شکل ٩‏ ابعاد هندسی صفحه‌ی فیلت تحت کشش یک طرفه 


ارائه آلکوریتم‌های بهینه‌سازی کارا با استفاده از روش رشد ... 


۳۹1 


وس 


ات ۲ 


شکل ۱۰ مدل اجزای محدود فیلت در نرم‌افزار 5 ۸5۷ 


۵ 4ووموو 
24 


7 


1000۸ 0۷ 


#4« «»مودد ۶ 


3033 ۹4 41 


شکل ۱۲ توزیع تنش ون مایزز در حالت بهینه 


جدول ۲تنش‌های ون مایزز در دو حالت بهینه نشده و بهینه برای فیلت 


نشریه مهندسی عمران فردوسی 


3 0۷۳۵۱ وععتاو عععع]-۷۵۵ 
رباص جع]۷ حیحصت صهع]۷ 
بهینه نشده 193 4عوو.0 096 
بهینه 134 2ع6د.! 068 


۹2 
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0 ۰1۰991 2-93 . 127 .-7 

7 6 1 . 7 1۰45 . 10۰07 . 6257 2-0196 
1007 294 > 44 وود 2304۰57 262 

(62 45 04 20 . 157 . (7 


(ب) 
شکل ۱۳ توزیع تنش اصلی اول الف) در اين پژومش, ب) مقاله[1۷۳۷]2 


در شکل (۱۳) توزیع تتش اصلی اول نسبت به 2 
نتایج ۱۷۷ [2]مقایسه شده است. نتایج نشان می‌دهد که 


در این پژوهش بعد از ۲ چرخه جواب همگرا 
می‌شود در حالی که در مقاله‌ی[۷۷۷]2 بعد از ۲۷ 


چرخه همگرایی حاصل شده است شکل (۱۶). ۱ ۰ ۱ ۲ 1 


0 5 10 15 20 25 30 


هم چنین» ماکزیمم تنش اصلی اول نرمال شده در این 5( 
مقاله ۱/۳۶ است که نسبت به نتایج [۷/۵]2 (۱/۵۲۱) 1 


۶ 260 0۵۲۵ ان ۱2۳ 
5 0۲۱۳۱۱۱ 


شکل ۳ سرعت همگرا شدن الگوریتم به‌کار برده شده برای 
درصد کم‌تر است. فیلت. الف) در این پژوهش, ب) الگوریتم [1۷۷]2 


سال بیست و پنجم. شماره دو. ۱۳۹۳ نشریه مهندسی عمران فردوسی 


۳۹ 


تحلیل اجزای محدود با مدول الاستیسیته ثابت 


تحلیل اجزای محدرد 


با مدول الاستیسیته تغییر بافته 


چرخه 


شکل ۱۵ ساختار الگوریتم حذف نرم 


روش رشد سازگار بیولوژیکی در بهینه‌سازی 
ساختار 


در بهینه‌سازی ساختار با استفاده از روش رشد سازگار 
بیولوژیکی از این رفتار درختان الهام گرفته می‌شود که 
در محل‌هایی که بار واردشده بیش‌تر است تنه‌ی 
درخت سفت‌تر و محکم‌تر است. در بهینه‌سازی یک 
سازه نیز می‌توان از این رفتار بهره گرفت به گونه‌ای که 
قسمت‌هایی از سازه که تتش بیش‌تری دارند از مواد با 
مدول الاستیسیته بالاتر و قسمت‌هایی که تنش کم‌تری 
دارند از مواد با مدول الاستیسیته کم‌تر انتحاب شوند. 
یکی از الگوریتم‌هایی که برای بهینه‌سازی ساختار 
سازه‌ها با استفاده از روش رشد سازگار بیولوژیکی 
ارائه شده. الگوریتم حذف نرم است. این الگوریتم بر 
پایه‌ی مدول الاستیسیته المان‌هاست و در هر مرحله 
مدول الاستیسیته المان‌های جسم را تابع تنش آن المان 
در مرحله‌ی قبل در نظر می‌گیرد. در نهایت المان‌هایی 
را که مدول الاستیسیته کم‌تری دارند و یابه عبارتی 
تنش کم‌تری دارناد را حذف و پا بی‌اثر می‌کند و 
الفان‌های سا منون السشست الا رامواار بسانه 
نگه می‌دارد. 
ساختار این الگوریتم همان‌گونه که در شکل (۱۵) نشان 
داده شده به صورت زير است. 
۱- محاسبه‌ی تنش‌های ون مایزز هر المان با در نظر 
گرفتن شرایط مرزی و بارگذاری و در نظر گرفتن 


نشریه مهندسی عمران فردوسی 


ارائه الکو ریتم‌های بهینه‌سازی کارا با استفاده از روش رشد ... 


مدول الاستیسیته ثابت برای همه‌ی المان‌ها. 
هر ول لامستیهه لمانها زا در تطر کتر وین 
تنش‌های ون مایزز به دست آمده در گام ۱ با 
استفاده از معادله‌ی (1). 
)0( (م.6 - ,150+ رط < ببوظ 


که رررظ و ,رت به ترتیب مدول الاستیسیته در مراحل 
1+ و 7 است و ,6 تنش ون مایزز در مرحله‌ی 
ام برای هر المان و مره تنش مرجع برای تمامی 
المان‌هاست. >افاکتور مقیاس است که در شروع فرآیند 
بهینه سازی می‌تواند ,6 / 100 انتخاب شود و در طی 

فرآیند افزایش یابد. 

۳- تحلیل سازه با المان‌هایی که مدول الاستیسیته آن‌ها 
تغییریافته است و محاسبه‌ی توزیع تنش. 

۶- تکرار مراحل ۲ و ۳ و حذف المان‌های بامدول 
الاستیسیته‌ی کم و در نهایت به دست آوردن 
ساختاری بهینه[14]. 

طراحی بهینه‌ی سازه‌ها تنها با استفاده از یک مدل 

اجزای محدود و بدون نیاز به تولید مدل جدید و یا 

مش‌بندی جدید به دست می‌آید. افزون بر این»عامل 
تعبین کننده در ساختار بهینه, تنش مرحم است و افزایش 
سیکل‌ها؛ که معمولا هزینه‌ی قابل توجهی به طراح 


نمی کند. 


کاربرد روش رشد سازگار بیولوژیکی در 
بهینه‌سازی ساختار 
مثال 1. در اين مثال یک تير یک‌سر گیردار مخروطی 
تحت بار عرضی در نظر گرفته شده است. در این مثال 
نیز از نرم‌افزار5 ۸۱5۷ برای تحلیل و پیاده‌سازی 
الگوریتم استفاده شده است. مدل اجزای محدود و 
شرایط مرزی و بارگذاری در شکل (۱۸) نشان داده 
شده است. ساختار بهینه‌ی تیر بادو تنش 
مرجع۱/۲(۳۵ و۳۷/۳۵ به دست می‌آید. همان‌گونه که 


سال بیست و پنجم» شماره دی ۱۳۹۳ 


در شکل (۱۷) دیده می‌شود. ساختارهای بهینه‌ی به 
دست آمده با تش‌های مرجع مختلف تفاوت 
چشم گیری دارند که نقش تعبینکننده‌ی تنش مرجع را 
نشان می‌دهد به گونه‌ای که با افزایش تنش مرجع» 
مقدار بیش‌تری از ماده حذف می‌شود و با کاهش آن 
حجم نهایی تیر افزایش می‌یابد. 

با مقایسه‌ی نتایج این پژوهش و نتایج 0 و 
همکاران[15] شکل (۱۸) مشخص می‌گردد که ساختار 
بهینه برای تیر مخروطی با تتش مرجع کم از نظر شکل 
مشابه ساختار به دست امده توسط ۹60 و همکاران 
است که با تغییرات ساختاری در فضای پارامتری و 
فننکین, یسو :تیار فسة کسمته: آ له اسست لسکا 
(۱۸- ج). با این تقاوت که حجم نرمال که عبارتست 
از نسبت حجم نهایی به حجم اولیه» برای سازه‌ی اینن 
پژوهش ۷ برای سازه‌ی ۹60 و همکاران ۰/۷ است 
شکل (۱۸- ب). از طرف دیگر ساختار بهینه برای تیر 
مخروطی با تتش مرجع بالاتر مشابه ساختار ٩60‏ و 
همکاران است که با بهینه‌سازی شکل تير مستطیلی به 


۳۷ 


دست آمده است شکل (۱۸-د). با این مزیت که حجم 
نرمال سازه‌ی این پژوهش در این حالت ۰/۳۸ و حجم 
نرمال ساختار 0 و همکاران است. 


مثال ۲. در این مثال یک تير با سطح مقطع متفاوت در 
نظر گرفته شده است. در شکل (۱۹) مدل اجزای 
محدود و شرایط مرزی و بارگذاری نشان داده شده 
است. با اعمال الگوریتم حذف نرم ساختار بهینه به 
دست می‌آید. همان‌گونه که در شکل (۲۰) دیده 
مي‌شبوده در اینجا نبز تشن مرجم عش تین کندده‌ای 
در ساختار بهینه دارد. در ضمن حجم نرمال برای تیر با 
ی اتکی مر 9 ی 
مرجع زیاد ۰/۲ است و این نتایج بسیار به‌تر از مثال 
قبلی و نتایج 0 و همکاران است. در نهایت می‌توان 
نتیجه گرفت که در حالت تنش مرجع بالا با شرط 
یکسان بودن شرایط مرزی و بارگذاری. ساختار بهینه‌ی 
تیرهای با شکل اولیه‌ی متفاوت. یکسان خواهند بود. 
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شکل ۱ مدل اجزای محدود تير مخروطی تحت بار عرضی 


تبث 


شکل ۱۷ ساختار بهینه برای تیر مخروطی تحت بار عرضیی الف) با تتش مرجع ۱/۲۵ ب) با تتش مرجع ۳۸/۵ 
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شکل ۲۰ ساختار بهینه برای تیر با سطح مقطع متفاوت تحت بار عرضی الف) با تتش مرجع ۱/۲/۸۵ ب) با تنش مرجع ۳/۳۵ 


الف 
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نتیجه گیری 

در این مقاله هدف به دست آوردن شکل و ساختار 
بهینه‌ی سازه‌های مکانیکی به منظور یک‌نواختی تنش» 
کاهش تنش‌های ماکزیمم و يا کاهش حجم است. با 
۱ 
شده است که شکل بهینه‌ی سوراخ و فیلت را به دست 
می‌آورد.در صفحه‌ی سوراخ دار. سوراخ دایره‌ای در 
شکل بهینه به صورت بیضی با قطر بزرگ‌تر در جهعت 
بار بیش‌تر تغییر یافته است. هم‌چنین» تنش متوسط در 
حالت بهینه ۲۲ درصد کاهش می‌یابد. سپس برای 
محاسبه‌ی عمر خستگی سازه با استفاده از باطاهکدی 
نوشته شده است تا عمر خستگی به صورت کانتور بر 
روی مدل نمایش داده شود. عمر خستگی در حالت 
بهینه تا ۶۰ برابر افزايش می‌یابد. در مورد فیلت نیز تنش 
متوسط در حالت بهینه ۳۰ درصد کاهش یافته است. 
میتی مافزیس شش تن آزن توسال هه در تن 
پژوهش نسبت به پژوهش‌های قبلی ٩‏ درصد کم‌تر شده 
که نشان‌دهنده‌ی مزیت روش پیشنهادی است. 

افزون بر این الگوریتم حذف نرم برای بهینه‌سازی 
ساختار تیر مخروطی و تير با سطح مقطع متفاوت اعمال 
شده است. نتایج نشان می‌دهد که با انتخاب تنش مرجع 
بالاتره ساختار بهینه برای تیرهای با شرایط مرزی و 
بارگذاری یکسان مشابه خواهد بود. اعمال این الگوریتم 
منجر به کاهش حجم نرمال و به دنبال آن کاهش وزن 
سازه نسبت به موارد مشابه پیشین شده است. 


فهرست نمادها و نشانه‌ها 
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